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Abstract 

The behaviour of the Dirac field in FRW space-time is investigated. 
The relevant equations are solved to determine the particle and energy 
distribution. The angular and radial parts are solved in terms of Jacobi 
polynomials. The time dependence of the massive field is solved in terms 
of known function only for the radiation filled flat space. WKB method 
is used for approximate solution in general FRW space. Of the two inde- 
pendent solutions, one is found to decay in time as the Universe expands, 
while the other solution grows. This could be the source of the local par- 
ticle current. The behaviour of the particle number and energy density 
are also investigated. It is found that the particles arrange themselves 
in a number and density distribution pattern that produces a constant 
Newtonian potential as required for the flat rotation curves of galaxies. 
Further, density contrast is found to grow with the expansion. 
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1 Introduction 

Since all matter is ultimately fermionic, behaviour of Dirac field in expanding 
Universe must reveal hitherto unknown features of the Universe. The properties 
of matter fields in general, and the massive Dirac fields in particular, must have 
profound consequences on the structure formation proecess. 

The spinor formalism developed by B. v. d. Waerdenpj in 1928 and later 
elaborated by O. Laporte, et al comprises all representations of the Lorentz 
group, even those not contained in ordinary tensor analysis. However, it was not 
widely used at that time. Later, in 1936, the spinor calculus was used by Diracjl] 
to derive wave equations for particles with spin greater than one half, that 
means wave equations for (elementary) particles that had not been discovered 
up to then. This made the spinor calculus known to a wider audience and Brill 
and Wheeler's analysis and application to the interaction of neutrinos with 
gravitational fields further added weights to this. When Newman and Penrose 
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developed a special formalismfS^ of projecting vectors, tensors and spinors onto 
a set of null tetrad bases (explained in detail in Ref. |Z])i it became a very useful 
tool to study the behaviour of quantum-fields of all spins in general relativistic 
back-ground gravity. This method has been used successfully in various black 
hole geometries jTJ [SJ [HI HJ and more recently to further study the behaviour 
of Dirac particles in different geometries |121 [T^ UM [T51 [T5] . 

In continuance of previous work |17L I18| on applying this method to inves- 
tigate the behaviour of the fields of electrodynamics in Friedmann-Robertson- 
Walker(FRW) space-time, we take up the Dirac field to get some insight on the 
process of structure formation. For our purpose, we first write the metric as 

[111120] 

ds^ = [dij"^ - dr^ - S"^ {dO^ + sin^ Odcf)'^)] (1) 

where a is the scale factor that depends on the conformal time rj which is 
related to the cosmic co-moving time thydt = adrj, and 

{sinr, K = 1 
K = . (2) 
sinhr, K = —1 

The flat case can be considered to be the limit of either the closed or open 
cases for if ^ 0. Thus, in the complexified r-plane, the closed Universe lies 
along the real axis, open along the imaginary and flat near the origin. So, we 
may consider the proposition that, as an overall flat space expands, it frag- 
ments into under-density voids that develop locally as an open FRW space 
along the imaginary axis, and the overdensity regions evolve locally into the 
large scale structures as the closed counterpart. We choose the null tetrad as 
If, = [1,-1, 0, 0], = ^^[1, 1, 0, 0], rrif, = -^[0, 0, 1, isine] and the com- 
plex conjugate to^so that the directional derivatives can be expressed as 



where. 



and 



^ d d . 

^ dr drj 
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The spinor equivalent of these directional derivatives are 
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dno' = D, dir = A, dov = S, and dw = S* 
Also, the non- vanishing spin coefficients [71 are given by (3 



(6) 
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2ii : 
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Dirac equations for spin-f/2 particles represented by a pair of spinors, P"* and , 



are 



and 



AB 



AB 



M - 



(7) 



(8) 



where cr*^^, are Pauli-niartices and M is the mass of the particle (expressed 
as the inverse of its Conipton wave-length). To write these equations in the 
Newman-Penrose formalism in a curved space-time, we replace the ordindary 
derivatives with covariant derivatives and Pauli-matrices by a generalized to 
curved spacetime as 



AB' 



1 

71 



Thus the Dirac equations are 



M 



£C'B' 



and 



where 



M -r' 
V2 



'AB'Q,^ + ^^P'^ec'B'^0, 



£C'B' 



1 

-1 



(9) 



(10) 



(11) 



(12) 



consider the case i?' = in Eq. (10 1: 



In order to write explicit forms of these equations in terms of spin coefficients, 

(13) 



0, 



Explicitly, in NP formalism this becomes 



dlO 



,P'')=i 



M 

71 



(14) 



Noting that = e'^^T^ = -Pq and P" = e°^Ti = Piand using ([g]), this 
equation can be written as 
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DP° + (Twoo'P'^ + Tnoo'P^) + S* - (rooio'^° + ToiicP^) - i^Q^' (15) 

Replacing the various spin coefficient by their special symbols (T], we obtain 

(I? + £-p)P" + ((5*+7^-a)Fl=^^Ql' (16) 
Similarly, for i?' = 1 in Eq. (|10[) gives another such equation and following 



the similar steps with the complex conjugate of Eq. (11) gives the other two 
equations. Hence the Dirac equations for the components of the Dirac spinors 
in NP formalism are given by 



(D + e- p)P'' + {6*+TT- a)P^ = i^Q^' , 

{A + fi- 7)pi + {S + P- r)P" = -^^Q°', 

(D + e* - p*)Q°' + iS + n* - a*)Q^' = 'i^P\ and 

V2 

- , - , M 

(A + ^* - 7*)Qi + [S* +13*- T*)Q'^ = i — P". (17) 

V2 

In the next section, the solution of these equations including angular, radial 
and time parts are presented. Section 3 deals with the particle current and 
the energy-momentum of the field. Concluding remarks are made in the last 
section. Although the solutions of Dirac equations have been obtained earlier 
|211l^[^[21] . they are either restricted to some special cases or the expressions 
are comparatively complex ones. With our choice of null tetrad bases, the 
separated equations are more tractable and the solutions can be written in 
terms of known mathematical functions. 



2 Solutions 



Substituting {V2SaP°, Sa^P\ -V2SaQ°\ Sa^Q^') with 

(0-i/2>"-i/2, 0i/2ia/2, 0-1/2^1/2, 01/2^^-1/2) in the Eqs. ([l7|, using the 
directional derivatives and spin-coefficients given in Sec. 1, the angular parts of 
Eqs. ( 17 1 are found to separate into 

1 



/:± + -cot(?)y±i/2 = TAr^i/2, 



while the radial and time parts satisfy 

r'D±(l)±i/2 = (/A T iMra) (j)^i/2 



(18) 



(19) 
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Substituting for I+1/2 = {C- + \ cot0) y_i/2from Eq. (18 1 into the one for 



Y_i/2, and a similar process with ^-1/2, gives us the two angular equations 
1 d . „ d . icos9 d'^ 1 2 n ^ \ 2 



■ ^7r^-7^o*^-o-^ P±i/2 = (20) 

^sm9o& oO sm 9 Oip^ 4 2 J ' 

Explicitly, the angular parts representing the two helicities are solved in 
terms of the Jacobi polynomials Pn"''^'' given by the spin weighted spherical 
harmonics 



±sYi"\e,cl)) = iVe^-^Cl - cos0)T±§(l + cos0)TTf p^("±«^™T«)(cos0) (21) 

of spin weight s = ^; here, the total angular momentum I — s + 1, s + 2, 
are half-integers as they already include the spin and the separation constant 
A = (I + 1/2). zbs^^/™ are appropriately normalized and weighted to give 
Ie=oId=o sYiV^ * sYJ^ sin 9 dddcj)— SuiSmm'- The radial and temporal parts can 



also be separated by writing R±T± ~ 4>i/2 i 4>-i/2 ii^ Eq. (19 1 to give the 
equations 

The separation constant k — pa can be identified with the comoving momentum. 
The solutions labeled ± can be decoupled easily [T7], and we believe that they 
represent the negative and positive energy solutions respectively. The decoupled 
equatins are of second order and are given by 

and 



dr S{r) J \dr ' S{r) 



^ ± iAIa^ ^ T iMa^ T± = ~k^T±. (24) 

2.1 Radial Solutions 

The resulting regular radial solutions 

i?± = A±(l-cosVKr] ' " (l + cos^/Kr] ' ' (25) 
P^^};^-^^^^ (cosVKr), 

are the appropriately weighted Jacobi polynomials that can be normalized to 
give \R±\^ dr = 1. The behaviour of R± is displayed in Fig. (jlj). It suffices 
to write the radial solutions for all K's in the form of Eq.(|25|, as the limit to 
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flat space is achieved with if 0, in which case the asymptotic behaviour of 
the Jacobi polynomial of large order leads to [25,] 



(26) 



where j„ are the spherical Bessel functions. It was shown in Ref. [T7] that 
R*^R- = —Rjf.R*_, so that A- = iA+. As the mass does not appear in the radial 
solution, the spatial solution of the massive and massless field equations are the 
same. The mass couples with the scale factor and appears only in the temporal 
part. One important aspect of the solution Eq 
has to be an integer for R± to be regular at r = and tt. As 
half-integer increasing in steps of unity from its minimum value, fc > | is also 
half-integer. Hence, the momentum of the free Dirac field in closed FRW space- 
time is quantized in steps of unity from | upwards. This minimum value of | 
can be interpreted as half unit of spin and one unit representing the mass in the 
lowest momentum state for the massive or one unit of lightlike momentum for 



(251 is that k - I -1 = 



> i IS a 
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the massless fields. Eq (22 1 leads to a number of Wronskian-like conditions, the 



most important of which that we will require later are 



IT+iVir_|' 



ik [r;r_ - t+ti] 

ik - \R\^ 



1 



(27) 



--j-RtR'- 

ar 



2.2 Time Evolution 

The decoupled temporal equation is 



drf 



k^ + M^a^ T iMa' T± = 0, 



(28) 



and a'2 = (^^^ = [^f^a'^ -I- (1 - f^A - ^ii - ^r) + ^Aia + ^r] is given 
by the Friedmann-Le Maitre equation with JIa, and fl^ being the densities 
contributed by the cosmological constant, matter and radiation respectively, in 
units of the critical density at the time tjq when the scale factor is normalized to 
oo = 1. Given a functional form for a (r?), it is always possible to write a series 
solution fT5] of Eq.(28l. It is instructive to staxt with solutions for the massless 
field. The two solutions are just 



— ikrj 



(29) 



where the multiplicative constants and their signs have been fixed by the con- 
ditions imposed by Eqs. ( p2| ) and (27 1. 

As shown in Ref. |17|. the time dependence can be solved in terms of a 
known special function (Whitaker) only in the case of the radiation filled fiat 
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Universe. This case is very important for the reason that it represents the nature 
of the early radiation dominated Universe from the very big-bang during which 
time, the density parameter was extremely fine-tuned to almost unity. So before 
trying to approximate the general solution, let us consider this in some further 
details. The solutions are better represented in terms of the parabolic cylinder 
functions 125] as 



T± iv) 



(30) 



where Eq.(22l has been used to relate the constants of the two independent 
solutions. \/ M Hq is the reciprocal of the geometric mean of the Hubble time of 
the Universe and the Compton time of the Dirac particle. From the discussions 
of the previous paragraph, we can fix the constants to be 



1 _ 
V2 



D 



, s k(l±i) 



_i(0) 



by requiring T± (0) ~ ±4^- The behaviour of the solutions Eq.(30l for some 
values of k are shown in Fig. O and p|) . 



Next, we try to approximate the solutions of Eq. (30 1 by WKB method 



assuming that a changes very slowly with rj so that a' (77) is small, i. e. the 
expansion rate of the Universe is small compared to the size of the Universe. 



Multiplying the derivative in Eq.(22l with S to incorporate the order of the 



approximation and substituting T± — exp 
get the equation for /: 



n=0 



E fndii) in Eq. (281 we 



(31) 



n=0 



^ m— 



Equating equal powers of 6 we get the equations for the respective orders of the 
approximations: 



/o 



v/fc2 + Af 



aeo, 



2/0/1 = -I (/o ± Ma') and 

n 

"*/n— 1 ~ ^ ^ /m/n— m- 



T?J. = 



(32) 



Solving to the first order, we find J fidi] — iln yl T f^, so we can write the 
solution to first order as 



T+ 



± 



2£n 



(33) 
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which resembles the solution in Minkowskian space-time. The solution also gives 
us another important information, viz., at late times (in the evolution of the 
Universe) when a — > oo, £o ^ so T+ — > and T„ — > —1. This shows that 
one of the solutions decays in time as the Universe expands, while the other 
solution grows. So, we suspect that the former one should be negative energy 
solution and the latter the positive energy solution because it is in consistent 
with the behaviors of the Dirac field in the Minkowski space where it is found 
that the negative energy solutions exist only in the relativistic regime while in 
the non-relativistic limit it is the positive energy solution that dominates |27| . 
However, it has to be confirmed with further investigations. 



3 Particle Current and Energy-Momentum 

The r] and r components of the four-current integrated over the angular co- 
ordinates are found to be (with S defined in Eq. Q) 



|i?^+|'|T+|V|i?-l'|r-l' 



= 47rS'''a'' + 71+) 

and Ana^rS^ = (TTT, - T+Tl) = |T+|t34) 

^ ' ik drj 

At this point, we are not much concerned with the other components, viz. c>c 
Y*Y_ - Y+Yl = and a'^J'^'S^ = -'-^^ {TIT_+T+T*_) {YIY_+Y+Y1) 
which is independent of (f). The equation of continuity becomes -^^na^S^ = 



— -§pCL^S^J^ — Tjr — %r-— ^^^^^3^ on using the relations Eq. (27 1. Integrating the 
comoving number density over the whole volume, the total number contained 
within the whole of the closed Universe gives us the normalization condition 
Nxa^ = dr J na^ sin^ rdn = |r+|^ + \T_f = 1 so that the total comov- 
ing number is conserved. In the extreme relativistic regime, i.e., in very early 
epoch of the Universe, when its size was small (a — >■ small), assuming a Fermi- 
Dirac(FD) type of distribution, and taking the volume of the closed space to 
he V — S^dr J dQ = 27r'^, we may write the average number density of the 
particles as 



= A E , = A £ ^^T^vf = 0.0827975. (35) 

27r2 ^ e'^+'+i + 1 27r2 ^ efe+3/2 i ^ ^ 

The FD value is na'^ = 0.0913454, and this distribution deviates significantly 
from Fermi-Dirac at low n. The number distribution, corresponding to the 
massless case, is shown and compared with the FD distribution in Fig. Q. 
On the other hand, integrating over the angular co-ordinates and only out 
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to a radius r, we have 



Na^ = I na^S^drdn 



\R-rdr 



R*,R- 



ik 



1"* order WKB 



Ma 



(36) 



for the relative number of particles within a finite volume. This WKB result 
are plotted for representative cases in Fig. ([5|. 

Particularly, for flat Universe, the comoving particle number and the particle 
current can be obtained using the expression for R± given by Eq. ( |26[) and 
T± given by Eq. ( 30 1 . Their sample plots have been shown in Figures 
(with different values of M) and ( 10 1 (with different values of conformal time rj). 



Next, we consider the energy-momentum tensor. At present we are con- 
cerned only with the 7777 and rr components giving the density p and the pres- 
sure Pr respectively. After integration over the angular co-ordinates, we end up 
with 



Ma{\Z^r-\Z\^ 



and 
AnS^ {p - 3P) 



k{\R.\' + \R\' 



{21 + 1) R*+R- 



smr 



(37) 



(r;T_ + T+T*) , 



-M[\Zj'-\Z\^ 



In the last equation for the trace, 3P = Pr + Pe + Ftp- Differentiating the energy 
density with respect to a we find the energy conservation law as 



47r S"—pa^ 
da 



(1^ 



-M |Z_ 



\Z\ 



AttS^ {p - 3P) a^ 



(38) 



WKB results for the quantities pa"* and Ea'^ — An S^pa'^dr are shown in Fig. 

(fTTl). 



Comparing Eq. ( 38 1 with Eq. ( 34 ) , it is important to note that 



AirS^na^ 



+ \R+T+\^ 



and AttS'^ 



d{pa*) 
d{aM) 



|i?_T_|^ - |P+r+|^ = AnS'^a^ (n_ 
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so that the rate of energy density change per unit mass is equal to the 
difference in the comoving particle number density of the '-' and '+' states. 

From Eq. (38 1, d [pa^) ~ —Pda^, which is the Second Law of Thermody- 
namics for adiabatic expansion. This equation can be processed to give the time 
independent entropy 

[p + P) 



S{r) ^ 
The total energy inside the closed Universe is 



(39) 



dr / pa sin rdfi 



^ -k{TlT^+T+T*_) - Ma[\T+\^ -\T^\^ 



(40) 



Dividing this by the conserved comoving total number Na^ — \, and using 



similar process with Eq. ( 38 1 , we get the equations for the average comoving 



energy of one Dirac particle and its time evolution 
ea = -k (T^T 



T+Tl] - Ma { iTi P - |T_n and 



dea 
da 



= -Ma [IT, \^ - IT \^ 



Eqs.(41 and 42) allow us to write 

2M\T±y 

k{TXT^ +T+T1] 



M =F -;— and 
da 

2 

da 



(41) 
(42) 

(43) 
(44) 



We can also express the evolution of the moving energy in any Friedmann type 
model by the equation 



a'^a'4- + 4fc2 + AM^aA ^ - Ul\ (ea) - 0. 
da da I da 



(45) 



Eq. (45) can be solved, at least numerically or by series method, for any be- 



haviour of a. 

For the flat Universe, the energy density and rate of flow of energy can be 
obtained explicitly as we did for particle number, by using the expression for 
R± given by Eq. (|26| and T± given by Eq. (ISOl) in Eq. |40| . The plots of these 



for some representative cases have been shown in Figs. (12), (13) and (14) 



1 diEa-^) d{Ea^) 



a di] 



da 



I {\T+\^ - l) r^ji-i/2ji+i/2 



M 
2k 



|A±|^^(l-2|r+|^) f {rj,+,/,Ydr 



(46) 
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4 Conclusions 



In this work, we have investigated the behaviour of the Dirac field in closed FRW 
spacetime using NP formalism - spinor analysis. This work has two important 
aspects. One is the method itself; we want to illustrate the fact that the NP 
formalism is a very powerful and convenient method to study the behaviour of 
spin fields in curved background spacetime. It yields more tractable equations 
for the eigenfuctions compared to the conventional tensorial method. It appears 
to be more natural as it goes into very quantal nature of the fields as illustrated 
by its ability to explain the absence of super-radiance of neutrinos in Kerr 
spacetime |7j. Another aspect is the solution of the field equations itself that 
have become more transparent in revealing the nature of the distributions of the 
particles and energy. Ultimately all matter is fermionic, so behaviour of Dirac 
field in expanding Universe must reveal the unknown features of the Universe 
and also unleash some facts about structure formation. 

While the free Maxwell and Dirac fields, in Friedmann-Robertson- Walker 
spacetime were investigated using the Newman-Penrose method, all the vari- 
ables were found to be separable, and the angular solutions were the spin- 
weighted spherical harmonics. The massless fields came to have the usual 
exponential time dependencies. All the radial parts reduced to the quantum 
mechanical barrier penetration problem, with well behaved potentials that are 
basically the centrifugal energies. The potentials seen by one component of the 
Dirac field, are interesting; its lowest angular momentum state sees no po- 
tential in the fiat universe, while it sees an attractive one throughout the open 
universe; from afar, all angular momentum states of this component see attrac- 
tive potential in the open universe. Consequences of this effect may provide a 
means to determine whether the Universe is fiat, open or closed. All the radial 
equations are solved. 

We have looked into the properties of the massive Dirac field in detail for 
closed FRW spacetime. The radial dependence in the open Universe is also 
the appropriate limit. The fiat case can be considered to be the limit of either 
the closed or open cases for if — >■ 0. Thus, in the complexified r-plane, the 
closed Universe lies along the real axis, open along the imaginary and fiat near 
the origin. So, we may consider the proposition that, as an overall fiat space 
expands, it fragments into under-density voids that develop locally as an open 
FRW space along the imaginary axis, and the overdensity regions evolve locally 
into the large scale structures as the closed counterpart. Hence the closed, open 
and fiat Universes represent different regions of a complexified space. 

The radial dependence is the same for both the massless and massive cases. 
The momentum is quantized in steps of unity from the minimum value of 3/2. 
The quantization condition for the comoving mementum of a gravitationally 
trapped Dirac particle is A; = / + n + 1, where n> is an integer, and the total 
angular momentum / is a half integer. 

The temporal solution (Eq. 33 I gives us an important information, viz,, at 
late times of the evolution / expansion of the Universe when a — > oo, Eq — > 
M, so T-)_ and T_ —1, showing that one of the solutions decays in 
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time as the Universe expands, while the other solution grows. This leads us 
to conclude that the former one may be the negative energy solution while the 
latter is the positive energy solution. This will be confirmed in future work on 
the representation of the Dirac matrices. The variations in number densities 
exist from the very initial time of the big-bang. There was indication from our 
early work |T7] that, although the comoving number density of the massless 
field is locally conserved at every spacetime point, it is not so for the massive 
case. Even though the Universe starts out with equal number of particles with 
|T_|_|^ = |T_|^ = 1/2, at late times the former goes down towards zero while 
the latter grows to unity keeping the sum conserved and equal to unity. As the 
universe expands and the particles lose their kinetic energy, the time dependence 
shows that the |Z+T+|^ representing the number of particles in this model will 
decay while |Z_T_ | grows. This may be mimicking the behavior of the massive 
Dirac field in Minkowski spacetime where the negative energy solutions exist 
only in the relativistic regime while the positive energy solution dominates in 
the non-relativistic regime. IIJ. This could be the source of the local particle 
current that exists at every point in the FRW spacetime. Nonetheless, the 
total comoving number within the whole volume is still conserved. Further, the 
currents are set up in such a way that the over-densities are further enhanced 
while the under-densities get more depleted as the Universe expands. 

The energy distribution also follows the number distribution. We find that 
the number and energy distributions that are setup have exactly the same form 
as required by the flat ratation curves of galaxies. The Newtonian gravitational 
potential that is required for the flat rotation curve is that the square of the 
rotational velocity of the particles, cx ^^^-^ oc ^^-^ = constant. Our results 

also show that the plots of the quantities ^^^^ and ^^^^ become flat at a short 
distance from the origin as the Figs. (|5]and 11 1 show. Even the surface densities 
of particles and energy at a distance r show the same nature. The particle 
number and energy distributions have also been plotted for different masses 
and different conformal times (in the units of Hubble time) and it is seen that 
density contrast is enhanced for larger mass M as well as with the elapse of time 
as expected. 

The over- and under-density fluctuations that exist from the very begining 
further get enhanced with the expansion of the Universe. As a consequence, 
fragmentation can be induced as the underdensity regions inbetween the over- 
densities become highly depleted. The angular-momentum contained by these 
structures must also be the source of rotation of galaxies and clusters. As our 
Universe is almost flat, any gravitationaly bound large scale structure will mimic 
a closed Universe. 

This work has to be continued further, and these preliminary results con- 
firmed by simulation. Also, quantization of the Dirac field must be done and 
the corrections accounted for conclusive arguments. Another extension would 
be the study of quantum electrodynamics of the interacting Dirac and Maxwell 
fields. 
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Figure 1: The radial dependence, R+ on the right and R^/i on the left. The 
I = 1/2, are filled for k = 3/2 and dotted for 5/2 respectively; the / = 7/2 are 
dashed for k = 11/2 and dot-dashed forfc = 17/2. 




Figure 2: The temporal dependence of the Dirac field in a radiation filled flat 
FRW space-time. r+ are shown by solid and T_ by dashed fines respectively. 
The upper two are real and imaginary parts for ^^^^ = |, and the lower two 

for ^ = 1 
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Figure 3: Solid curves are|T+|^, and dashed are|r_|^; ^^^^ = | are at the left 

and ^^^^ = i at the right. \T^f are seen to decay to zero while|T_|^ grow to 
unity. 
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Figure 4: The relative particle number distribution of extremely relativistic 
fermions in closed FRW space. The momentum k is quantized upwards in steps 
of unity from the minimum value of 3/2, and the filled dots represent the relative 
number of particles with that momentum. The filled curve is the Fermi-Dirac 
distribution. We see that at low momenta, the distribution is quite different, 
but for high k the two start to coincide. 
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Figure 5: The number of Dirac particles in closed FRW space-time. The upper 
two are the local surface densities of particle number, and the lower two are 
the total number divided by r. On the left are k = 3/2, I = 1/2, and on the 
right are k = 17/2, / — 7/2. Solid curves are the respective quantities at the 
instance of the big bang when a = 0, and the dashed at a = 10. It is seen 
that as the Universe expands, the density contrast become more enhanced. The 
lower graphs show that the particle distribution has the same behaviour as that 
demanded by the flat rotation curves of galaxies. 
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Figure 6: The comoving particle number (Na^/r) for flat Universe at 77 — > 0. 
The normal solid line, dashed line and thick solid line are for k = 3/2, 5/2 and 
11/2 respectively, all of which are for I = 1/2. They arc appropriately scaled for 
comparision purpose. The quantities are found to saturate to a constant value 
resembling to the flattening of the rotation curves of galaxies. 
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Figure 7: The comoving particle number (Na^/r) for massive flat FRW space- 
time with I = 1/2 and k = 3/2. The normal solid, dashed and thick solid curves 
are for M =0.01, 1 and 100 respectively. The behaviour is same as in the case of 
masslcss field except that amplitude of fluctuations are larger, and increase for 
increasing mass. The graphs for the three cases have been appropriately scaled 
to make of almost equal heights for comparision. 
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Figure 8: The comoving particle current for given mass and conformal time in 
flat FRW space time. The thick soHd curve is for i = 1/2 and k = 3/2; the 
dashed curve is for I = 3/2 and k = 7/2; and normal solid curve is for I = 
5/2 and k = 11 2. The lower one is the particle number divided by r. The 
later graphs show that the particle distribution has the same behaviour as that 
demanded by the flat rotation curves of galaxies. 
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Figure 9: The rate of flow of comoving particles through a unit surface at r in 
flat FRW space time. The thick solid curve is for / = 1/2 and k = 3/2; the 
dashed curve is for / = 3/2 and k = 7/2; and normal solid curve is for I = 5/2 
and k = 11/2. 
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Figure 10: Evolution of Na^/r at different conformal times 77 for M = 1. The 
normal solid curve is for ^/MHqT] = 1, the dashed curve is for ^/MHorj = .01 
and the thick solid curve is for y/MUQ-q = 10, 000. All the graphs are for I = 
1/2 and k — 3/2. It is seen that the density contrast increases with the elapse 
of time. 
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Figure 11: The energy of the Dirac field for the same values of k and I as in Fig. 
([5]). But here the values are scaled for a = and a = 10 to appear to be almost 
equal so that they can be compared. The energy follows the particle number as 
it should. 
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Figure 12: Comoving energy flow rate for flat universe at a given time and 
mass. Thick solid, dashed and normal solid lines are for k =3/2, 5/2 and 7/2 
respectively. All are for I = 1/2. 
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Figure 13: Rate of flow of comoving energy for the flat universe. The upper 
one is the energy divided by r and lower one is the energy flowing out of a unit 
surface area at a distance r. All are for 1 = 1/2 with k =3/2, 5/2 and 7/2 
represented by thick solid, dashed and normal solid lines respectively 
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